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ABSTRACT 








It is known that the flowfield of the SRM (Solid Rocket Motor) is very complicated due to the 
complex characteristics of turbulent multi-phase flow, chemical reaction, particle combustion, evaporation, 
breakup and agglomeration etc It requires multi-phase calculations, chemical reaction simulation, and particle 
combustion, evaporation, and breakup models to obtain a better understanding of thermophysics for the SRM 
design using numerical methods. Also, the slag buildup due to the molten particles is another factor affecting 
the performance of the SRM Thus, a more realistic simulation is needed to provide a better design guide to 
improve the performance of SRM To achieve this goal, the VOF (Volume Of Fluid) method is used to 
capture the free surface motion so as to simulate the accumulation of the molten particles (slag) of SRM A 
Finite-rate chemistry model is used to simulate the chemical reaction effects. For multi-phase calculations, 
Hermsen combustion model is used for the AL particle combustion analysis and Taylor Analogy Breakup 
(TAB) model is used for the particle breakup analysis An interphase mass-exchange model introduced by 
Spalding is used for the evaporation calculation The particle trajectories are calculated using a one-step 
implicit method for several groups of particle sizes by which the drag forces and heat fluxes are then coupled 
with the gas phase equations 

The preliminary results predicted a reasonable physical simulation of the particle effects using a 
simple 2-D solid rocket motor configuration It shows that the AL/AL203 particle sizes are reduced due to 
the combustion, evaporation, and breakup. The flowfield is disturbed by the particles Mach number 
distributions in the nozzle are deformed due to the effect of particle concentrations away from the center line 

The RSRM (Redesigned Solid Rocket Motor) geometry at 67 seconds is employed to investigate the 
slag behavior in the aft-end cavity with the combustion, evaporation, and breakup models. The particulate 
phase was assumed to be aluminum oxide (AL203) for the preliminary study. It is assumed that the 
propellant grain of the aft-end cavity has burned out completely at 67 seconds. The geometry and mass flow 
rate information were provided by the NASA Marshall Space Flight Center. The slag may flow out of the 
cavity and enter the nozzle due to the accelerations The molten particles entering the aft-end cavity merge 
with the slag. The volume of the slag will grow and affect the performance of the RSRM. This shows that the 
effects of particles and slag on the flowfield are very significant From the caculation, a flow vortex exists in 
the aft-end cavity of the RSRM. A stagnation point on the wall is captured. This flow impingement may 
cause the erosion of the wall. The shape of the vortex is changed due to the slag. The particles entering the 
cavity may become slag and either flow' into or out of the cavity depending on the temperature and the surface 
tension of the molten particles. An axial gravity force of 2.4g is assumed to simulate the RSRM flowfield at 
67 seconds. 

The flowfield analysis using the FDNS code in the present research using the proposed models should 
provide a design guide for the solid rocket motors The obtained results can give the designer a basic guide 
line for the use of materials and the nozzle geometry to improve the performance of SRMs 


1293 



Numerical Investigation of Slag Behavior for RSRM 




T3 

C 

59 a 

bD | — 1 

_ (*> 
to qi 

A o 
(Z> S3 
. 



as 

S3 

a> 

JS 

V 

<z5 

> n 

£ 

3 

a-* 


OD 

S3 

• M 

La 

a> 

S3 

‘3d 

S3 

w 


'O 

e 

C3 


La 

<u 

S3 

0> 


U 



a > <d 
+* o 

«) S 

• Q- 
W (Z) 



o 

O 

o 



n 

fO 

O 

w 


1294 


13th Workshop for CFD Applications in Rocket Propulsion 



OBJECTIVE 
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VOF Model 

The VOF transport equation is given below: 


da 

~a 



da 

ck { 



where a = 1 stands for liquid and a = 0 is for gas. 
The interface is located at 1 > a > 0. For a given 
solution of a field, equation (6) can be recast as: 


for compressible gas : 

d Pm <t> &*(«-«'),* 
dt + dx i 

for incompressible gas : 
d<j> 


= S^a < 0.01 


dt 

and 


+ P 


m 


(«-««), 


d<k 

-f=S<,a > 0.0 


p m = Max\p g ,ap 


The interface a solution compression procedure is 
expressed as: 

1 > 0.5 + f(a old - Of)]} 


new 


= Max 


0 , Min 
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and 



[interface volume ) 


new 


[interface volume) 


initial 


The surface tension forces in the continuum 
surface force model is formulated as continuous 
body forces across the interface. These forces can 
be written as: 


A 


F. = - <j V n 
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F.= - <j\ V n 
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jCty + 
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v y j 


,for2D,axisymme trie case only 


A 


F z = - cr|^V n Ja x , for 3 D case only 
w here 

= surface tension constant 


A A A A 

Vn = a xx + a yy + a 


a is 0.5 for the free surface. The VOF method is 
used to represent the tracking of the free surface 
between the liquid and gas phase. 
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APPLICATIONS 
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2D Grid system of RSRM configuration. 


XMIN= 1.2SE+02 
XMAX= 1 5S1-H02 
YMIN=— 8 40E-K2)0 
YMAX= 1 42E+01 



Mach number contours of RSRM, no particles. 


Co 1 o i — ''lap 
a 0 0000E400 
b 0 3244E-02. 
c 1 6S4GE-01 
d 2 4973E-01 
e 3 32S7E-01 
f 4 . 1622E-01 
g 4 994EE-01 
h 5 0270E-01 
i S 05939-01 
j 7 4319E-01 
k 0 3244E-01 
1 9 1560E-01 

m 9 9GS3E-01 
n 1 2621E-H00 
o 1 1054E+30 
p 1 248EE-K20 
q 1 3319E+00 
r 1 4151E-H00 
s 1 4903E+00 
t 1 501GE-H30 
u 1 S04€E+00 
v 1 7461E-H20 
w 1 S313E-H20 
x 1 9140E+00 
y 1 997GE+00 
z 2 0811E+00 


1301 




T 



Slag accumulation in the aft-end cavity 



91531 sec. T = 2.31444 sec 






Slag accumulation at T = 2.3144 sec 
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Slag Accumulation for 3D RSRM at T = 0.031687 
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Slag Accumulation for 3D RSRM at T = 0.031687 
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Slag Accumulation for 3D RSRM at T = 0.031687 
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Slag Accumulation for 3D RSRM at T = 0.031687 



1309 


Slag Accumulation for 3D RSRM at T = 0.031687 
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